MHD Turbulence Simulation in a Cosmic Structure 

Context 

T.W. Jones^'^, David H. Porter^, Dongsu Ryu"^ and Jungyeon Cho^ 

^School of Physics and Astronomy, 
University of Minnesota, Minneapolis, MN 

^Minnesota Supercomputing Institute, 
University of Minnesota, Minneapolis, MN 
^Department of Astronomy and Space Science, 
Chungnam National University, Daejun, Korea 

August 8, 2011 
Abstract 

The gaseous media of galaxy clusters and cosmic filaments, which 
constitute most of the baryonic matter in the universe, is highly dy- 
namic. It is also probably turbulent, although the turbulence prop- 
erties are poorly known. The gas is highly rarefied, essentially fully 
ionized plasma. Observational evidence suggests intracluster media 
(ICMs) are magnetized at some level. There are several possible ori- 
gins for ICM seed fields; the observed fields are likely the result of 
turbulence in the ICM. We are engaged in a simulation study designed 
to understand in this context how very weak initial magnetic fields 
evolve in driven turbulence. We find that the magnetic fields even- 
tually evolve towards equipartition levels with the vortical, solenoidal 
kinetic energy in the turbulence. As they do so the topology of the field 
structures transition from filamentary forms into ribbon-like structures 
in which the field orientations are laminated with vorticity structures. 



1 Introduction 

Most of the baryonic matter in the universe is outside stars and galaxies. It 
exists primarily as very diffuse plasma in large scale cosmic filaments and 
galaxy clusters that have formed and are still forming through gravitational 
collapse of primordial density fluctuations. Both theory and observation 
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have established that the diffuse intra-cluster media (ICMs) are highly dy- 
namical environments with active "weather" driven by ongoing accretion 
and merging activity, as well as large energy inputs from starburst-driven 
galactic winds and very fast outflows from active galaxies. The ICMs are 
criss-crossed by complex winds at a fair fraction of the local sound speed that 
generate weak-to-moderate-strength shocks, contact discontinuities (known 
in the community as "cold fronts" ) and bulk shear. Provided Reynold num- 
bers are large enough such flows should become turbulent. Turbulence is 
clearly manifested in simulations of cluster formation, and there is growing 
observational support as well. For instance random flow velocities appear to 
have reduced resonance scattering of the 6.7 keV Fe emission line [1] of the 
ICM in the Perseus cluster. Thermal pressure fluctuations in the Coma clus- 
ter are consistent with Kolmogorov turbulence [2] . Patchy Faraday rotation 
patterns looking through several clusters also indicate highly disordered and 
possibly turbulent magnetic field structures [3] , as does the absence of large 
scale polarization in the diffuse synchrotron emission of radio halos seen in 
a growing number of cluster cores (e.g., [5). Direct information about the 
existence of turbulence in cosmic filaments is currently lacking, although 
there are theoretical reasons to expect the diffuse media of filaments also to 
develop turbulence (e.g., [5]). 

Turbulence in these environments is important to understand for many 
reasons. Turbulent pressure provides support against gravity, so is rele- 
vant to cluster mass estimates made from X-ray measurements. Turbulence 
transports entropy, metals and cosmic rays; all important cluster evolution 
diagnostics. It transports and amplifies magnetic fields (focus of this paper), 
which in turn control the viscosity, resistivity and thermal conduction in the 
diffuse media, as well as the propagation and acceleration of cosmic rays. 
There is, of course, an extensive literature on turbulence, including MHD 
turbulence (see, e.g., [6] for a review). Much of the astrophysical MHD tur- 
bulence literature (e.g., [7]) aims to understand galactic, interstellar media, 
where the magnetic fields are relatively stronger (/3 = Pg/Ps ~ 1) than 
in cluster and filament media (/3 >> 1). The simulation effort behind this 
report focuses on driven MHD turbulence in situations where the initial 
magnetic field is very weak, so especially relevant to clusters and filaments 
(see [1] for an early report from this study). 
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2 Generation of Turbulence and Magnetic Fields 



2.1 Turbulence and Vorticity 

As noted above the environments of interest are commonly filled with strong 
drivers of gas motion. Turbulence describes motions possessing significant 
random velocities. This random velocity field can include both compressive 
(V • n 7^ 0) and vortical, or solenoidal (oj = V x -u 7^ 0) components, 
where uj is vorticity. Significant amplification of magnetic fields depends 
on the presence of fiow stretching, so on the vortical velocity component 
(see equation E]). Thus, understanding MHD turbulence begins with an 
identification of the sources of vorticity and the manner in which vorticity 
evolves. 

The equation of motion for a viscous fluid can be expressed in terms of 
vorticity as, 

^ = V X (n X (j) + i^V'^uj + -^Vp X VP, (1) 
at 

where v is the kinetic viscosity (assumed constant and isotropic) [9j. In an 
ideal {y = 0) flow in which the baroclinic term, Vp x VP, vanishes the 
net vorticity of a fluid element is conserved {d/dt J uj ■ da = 0). When the 
pressure and density gradients are not aligned, such as in colliding flows, 
vorticity can be added through this term, while viscosity leads to decay 
of vorticity. In truly isothermal flows, such as those in the simulations we 
report here, the baroclinic term will always be zero, since P cx p everywhere. 

On the other hand a steady, uniform fiow (w = 0) obliquely crossing a 
curved shock surface will exit downstream with a post-shock vorticity given 
by "Crocco's Theorem" |10j . 



wcs = ^-^^^^KtJi X h, (2) 
P1P2 

with pi and p2 the upstream and downstream gas densities, Ui the upstream 
fiow velocity in the shock rest frame, K the shock surface curvature tensor, 
and h the shock normal unit vector. Crocco's theorem depends only on 
mass and momentum conservation at the shock (i.e., on differential fiow 
refraction), so it applies even in isothermal shocks, where baroclinicity is 
absent. We will explore this vorticity source explicitly below in simulations 
of compressible, isothermal turbulence driven entirely by sound waves. 

Turbulence develops as motions driven on a scale, L^, cascade into 
chaotic motions on smaller scales, provided the viscous dissipation scale. 
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Ivisc, is much smaller than L^- In our context the driving scale is generally 
comparable to such things as the curvature radius of a shock, the size of a 
cluster substructure core, or the scale of an active galaxy or starburst wind 
outflow. These likely range for clusters from ~ 10s of kpc to ~ 100s of kpc. 
In filaments even larger driving scales are likely. The appropriate viscous 
dissipation scales in these media are far less clear. They are hot, ionized and 
very diffuse, so Coulomb collisions are ineffective. The associated mean free 
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path, Xcoui ~ 1 kpcTj,gv/('^-3^th,ioo)5 ranges from 10s of pc to 10s of kpc in 
these environments. Here, T^eV is the plasma temperature in keV, n_3 is the 
density in 10~^cm~^, and uth,ioo is the ion thermal velocity in 100 km/sec. 
The corresponding viscosity, ~ uth^Coui is very large, and associated the 
Reynolds numbers. Re ~ LaU/v ~ few x 10, rather small with U ~ uth the 
flow velocity on the driving scale. Hence the viscous dissipation scale due 
to Coulomb scattering alone, Imsc ~ Ld/J^J^ , would range from fractions of 
a kpc in cool cluster cores to at least several 10s of kpc in cluster outskirts. 

On the other hand, the presence of even a weak magnetic field may 
reduce the dissipation scale as a result of gyroscale instabilities, such as 
the firehose and mirror instabilities. Then the scattering of particles by the 
resulting magnetic fluctuations could reduce the particle mean free paths, so 
also the viscous dissipation scale [ITj. The detailed picture is still uncertain, 
however. We assume below that the physical dissipation scale is at least as 
small as the effective dissipation scale of our "ideal fluid" simulations; i.e., of 
order the grid resolution. In cluster contexts this would correspond roughly 
to kpc dissipation scales. 

2.2 Magnetic Field 

The magnetic field evolution in a conducting medium is governed by the 
induction equation. For a generalized Ohm's law in the MHD approximation 
this is (e.g., [laila]) 

dB 1 

_ = V X (nx ^) + r/V2^ ^Vne x VP^, (3) 

at en^ 

where rj is the resistivity (assumed constant and isotropic,) while rig and Pe 
are the electron density and pressure, respectively. The last term in equation 
[3l which is analogous to the baroclinic source term for vorticity, comes from 
different electron and ion mobilities. It is the so-called "Biermann Battery" 
source term for magnetic fields frequently invoked as a contributor to the 
generation of cosmic magnetic fields, especially at curved shocks. It is not 
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commonly included explicitly in MHD simulations and not in those discussed 
here (but see, e.g., [5l [T2]) 

The first term on the right of equation [3] accounts for the important 
"stretch-fold" mechanism that leads to "small scale" or "turbulent" dy- 
namo field amplification. The resistance term controls field dissipation, of 
course. The resistive dissipation scales. Ires, are also uncertain in these envi- 

1/2 

ronments. In a turbulent flow with r] « v, we would have Ires ~ Ld/Rm , 
where Rm ~ LdU /rj is the magnetic Reynolds number. It is possible that 
the magnetic Prandtl number, i-*r,m = Rm/Re = v/r] > 1 in the media 
of interest here. To illustrate, if one assumes Coulomb scattering controls 
both dissipative processes, Pr,m >> 1 (e.g., [H]). In the simulations re- 
ported below the resistivity, like the viscosity has numerical origins, so that 
the resistive dissipation scale is also similar to the grid resolution; thus, 

Pr,m ~ 1- 

3 Simulation of Cosmic-Scale MHD Turbulence 

Isolating turbulence from coherent "weather" in the very complex and inho- 
mogeneous flows associated with clusters and filaments is difficult (but see, 
e.g., O IH]). As an alternative tool to explore some of the basic physics we 
are conducting a high resolution simulation study of the evolution and satu- 
ration of driven 3D MHD turbulence in computational domains that resem- 
ble these media. Since these media, while magnetized, are not magnetically 
dominated, we focus on turbulence developed with initially very weak mag- 
netic fields. The full study considers both compressible and incompressible 
fluids as well as ideal and non-ideal media with a range of magnetic Prandtl 
numbers. We discuss here, however, only some cases of ideal, compressible 
flows in isothermal media. The simulations used an isothermal ideal MHD 
code that is an updated version of the code presented in [15]. Initially the 
medium has a uniform density, p = 1, gas pressure, Pg = 1 (so isothermal 
sound speed, Cg = 1) and a uniform magnetic field with /3 = Pg/Ps = 10^. 

The cubic box has dimensions Lq = Lx = Ly = = 10 with periodic 
boundaries. The box sound crossing time is thus 10 units. Turbulence is 
driven by velocity forcing drawn from a Gaussian random field determined 
with a power spectrum, P^ oc ex.p{—8k / kpeak) , where kpeak = 2/co (^o = 
27r/Lo), and added at every At = O.OlL/cg. The power spectrum peaks 
around ka ?a 1.5ko, or around a scale. La 2/3Lq. The amplitude of the 
perturbations is tuned so that urms ~ 0.5 or Mg = urms/cs ~ 0.5 at 
saturation, close to what results in full cluster simulations (e.g., [5| \T6\ [T9]. 
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Figure 1: Left: Evolution of kinetic, Ek, and magnetic, Eb, energies in 
simulations of ideal, compressible solenoidally driven (/^ = 1) MHD turbu- 
lence for two grid resolutions. Right: Power spectra, E{k), of kinetic and 
magnetic energies at t = 20 and t = 130 in the 2048^ zone simulation. 

Using a Helmholtz decomposition, the driving velocity field is separated 
into solenoidal (V • 5m = 0) and compressive (V x Su = 0) components. The 
fraction of the total driving kinetic energy put into solenoidal motions is 
designated below by the symbol, fs- 

Results are presented here for purely solenoidal, = 1, and, for com- 
parison, purely compressive, /s = driving. Turbulence driving in a cosmic 
structure formation context will be somewhere between these two extremes. 
They are useful, however, in allowing us to focus on the the connection be- 
tween magnetic field amplification and vorticity on the one hand, and on 
the possible roles for shocks on the other hand. Our preliminary analysis 
of intermediate cases, < < 1, indicate behaviors that one might rea- 
sonably intuit from the two extremes. For the /s = 1 case we show results 
from simulations carried out on both 1024^ and 2048^ grids. While not fully 
converged they agree well in their general properties. For compressive driv- 
ing, fs = 0, we present results from a 512^ simulation, which is our highest 
resolution run to date for that case. The turbulent kinetic energy balance 
seems relatively well converged at this resolution , as well as the kinetic 
energy power spectra on large scales. The magnetic field properties may not 
yet be converged for reasons outlined below. 

Our setup gives a characteristic timescale of bulk motion, td = L^/urms ~ 
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Figure 2: Magnetic energy density distributions in solenoidally driven MHD 
turbulence. Left: Log{EB) at t = 20 in the 2048^ simulation. Right: 
Log{EB) at t = 130. "Cool" is weak; "hot" is strong. 

15. In that time the motions in both cases spawn some form of hydrody- 
namical turbulence with power from the driving scale down to the viscous 
dissipation scale (a few grid zones). Consider first the case with solenoidal 
driving, which is illustrated in Fig. [TJ It shows that the mean turbulent 
kinetic energy density, Ek grows and peaks at t ~ t^, with a value corre- 
sponding to urms ~ 2/3. Subsequently, Ek slowly declines as the mean 
magnetic energy density, Eb, grows, the kinetic energy power spectrum, 
Exik), at t = 20, also shown in Fig. [H exhibits a peak at k/ko ~ 2, near 
the driving scale. It takes a Kolmogorov-like, inertial form, E^ik) oc k~^l^ 
for k/ko < 50. By this time energy has cascaded from the driving scales 
far enough that the motions, with still-negligible magnetic backreaction, are 
reasonably described as classical, hydrodynamical turbulence over a modest 
range of scales. The kinetic energy is predominantly solenoidal; the ratio 
of solenoidal to compressive kinetic energies at saturation, Ek,s/Ek,c ~ 15. 
Consequently, the compressive motions play almost no role in this case. In- 
deed, the properties of analogous incompressible turbulence simulations are 
very similar. 

Once turbulence develops, both vorticity and magnetic energies undergo 
inverse cascades from small to large scales, with the coherence lengths of 
their filaments growing accordingly. This is evident for the magnetic field in 
the power spectrum changes in Fig. [TJ The inverse cascade of magnetic en- 
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ergy can be understood as follows. The field is wrapped more quickly around 
smaller scale eddies, because the eddy turn over time varies as ti oc P^^. 
Maxwell stresses, oc (V x B) x B, then, feed back on the kinetic turbu- 
lence, causing significant modifications in the fluid motions, thus saturating 
the magnetic fleld growth on a given scale, I, when Eb{1) ~ Ex{l)- Since 
the turbulent kinetic energy on a scale Ek{1) oc the saturation scale 
of the magnetic turbulence should evolve over time as oc t^^^, while the 
magnetic energy grows as Eb oc t, both consistent with Fig. [TJ Eventually, 
as Ib approaches the outer scale of the kinetic turbulence, L^, the scalings 
break down and turbulence reaches saturation where the ratio of the total 
magnetic to kinetic energy is Eb/Ek ~ 2/3 (see also, e.g., [20l EH [22]). 
Neither the kinetic nor the magnetic energy power spectra, nor their sum 
are well described as Kolmogorov in this saturation state. 

We also emphasize an interesting topological transformation in the flow 
structure as turbulence proceeds through the linear growth to the saturation 
stage. Fig. [2] displays the different topologies of the magnetic flux structures 
at t = 20 and t = 130. At the earlier time the fleld is organized into individ- 
ual fllaments. At the later time those fllaments have evolved into striated, 
ribbon-like forms (see also [21]). Close examination reveals the ribbons to be 
laminated, with magnetic fleld and vorticity interleaved through each cross 
section on scales of the order the dissipation length. In hydrodynamical 
turbulence such ribbons would be unstable, but the interplay of vorticity 
and magnetic fleld seems to stabilize them in MHD turbulence. 

The compressively driven turbulence, fs = 0, develops rather different 
properties, as illustrated in Fig. [3l The plots show solenoidal, compressive 
and magnetic power spectra in this case at t = 625 (right side), and for 
comparison the same power spectra in the previously discussed, /s = 1 case 
at t = 130 (left side). The /s = kinetic energy spectra resemble a Burgers 
scaling, Exik) oc k~'^, rather than a Kolmogorov scaling. That steeper 
scaling results from the dominance of shocks in the turbulence (e.g., |23)). 
despite the fact that the rms velocities in this turbulent flow are subsonic. 
It is not surprising then that the turbulent kinetic energy is predominantly 
in compressive modes, with the energy ratio, Ek,s/Ek,c ~ 1/15, the inverse 
of our /s = 1 result. 

Given the absence of any vorticity in the driving and the absence of any 
baroclinic vorticity sources, it is actually remarkable that there are vorti- 
cal motions at all. As noted earlier, the vorticity in this case is created at 
curved shocks and especially intersecting shocks in accordance with Crocco's 
theorem (equation [2]) . This effect can be seen clearly in the simulation, es- 
pecially early on, when the motions are essentially all compressive. Fig. 
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Figure 3: Energy power spectra in isothermal MHD turbulence with 
solenoidal driving (/^ = 1, left) and compressive driving (/^ = 0). Ki- 
netic energy in solenoidal motions, EK,s{k), compressive motions, EK,cik), 
their sum, EK,t{k), along with magnetic energy, Esik), are shown. 

m shows a small, planar slice from this simulation at t = 5, just as inter- 
secting shock structures are first forming. The associations between shock 
structures, vorticity and relatively stronger magnetic fields is obvious. 

Since the magnetic field amplification depends on the vortical motions, 
which are an order of magnitude smaller in the /s = case at hydrody- 
namical saturation than in the /s = 1 case, it is not surprising that the 
magnetic field grows much more slowly and is much weaker in the compres- 
sively driven case. In the /s = 1 case the linear growth of the magnetic field 
saturations around t ~ 80 — 100. In contrast, for the fs = 0, compressive 
driving case, the magnetic field is still in the linear growth phase at the 
end of the simulation, t = 625. The ratio of magnetic to solenoidal energies 
is Eb/Ek,s ~ 1/4, so still shy of saturation. Just as for the fs = I case, 
saturation of magnetic energy should develop when Eb/Ek,s ^ 2/3. This 
delayed and reduced growth is made more exaggerated in comparison to the 
steeper kinetic energy power spectrum, since small scale eddy turnover times 
(t/ oc l/u{l)) are longer with Burgers scaling. This also makes the magnetic 
field generation more sensitive to the grid resolution, since it depends on an 
inverse cascade, and so depends on solenoidal power on small scales. 
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Figure 4: Slice at t = 5 in the fs = simulation, revealing generation 
of vorticity and associated amplification of magnetic fields at intersecting, 
curved shock surfaces. Left: V • u, identifying shock structures. Center: 
magnitude of vorticity, \uj\. Right: Magnetic field strength. 

4 Conclusion 

Processes such as shocks and high speed outflows are likely to drive tur- 
bulence in the diffuse media in galaxy clusters and cosmic filaments. The 
detailed physics is difficult to model analytically. However, simulations al- 
low us to explore it in some detail. Magnetic fields are integral components, 
both in the microphysics of the media and in large scale dynamics, even 
though the fields themselves are not likely to be dynamically dominant. 
The magnetic fields also potentially provide critical diagnostics of the media 
and their dynamical states. The development of the magnetic fields through 
turbulence is likely, although the details of that development depend again 
on the microphysics (which itself is controlled by the magnetic fields) and 
the processes that drive the turbulent motions. Growth of an initially very 
weak field to saturated balance with vortical turbulent motions requires 
timescales long compared to the hydrodynamical turbulent timescale. 
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